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Abstract 

We consider the ferromagnetic Ising model with Glauber spin flip 
dynamics in one dimension. The external magnetic field vanishes and 
the couplings are i.i.d. random variables. If their distribution has 
compact support, the disorder averaged spin auto-correlation function 
has an exponential decay in time. We prove that, if the couplings are 
unbounded, the decay switches to either a power law or a stretched 
exponential, in general. 

1 Introduction and main results 

In one dimension the Ising model with spin flip dynamics has exponentially 
fast mixing in time, as is reflected by the fact that the self-adjoint generator 
of the stochastic dynamics has a spectral gap, see ^ for example. One 
might wonder what happens to the exponential decay when the couplings 
are disordered. If the couplings are uniformly bounded, it is proved in 0, |^ 
that the generator still has a spectral gap. Thus the case of interest is when 
the couplings are unbounded. It is easy to see that then the spectral gap 
vanishes with probability one. The goal of our paper is to estimate how the 
missing spectral gap is reflected in the decay of the disorder averaged spin- 
spin correlation. In particular we will have to identify those realizations of 
the couplings which are responsible for a slow decay. 
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The model under the study is the one- dimensional Ising model with formal 
Hamiltonian 



^xO'x-lO'x- 



(1) 



Here = ±1 are the Ising spins, a E Q = {1,-1}^, and uj^ are the couplings. 
We assume that {a;^., x G Z} are i.i.d. random variables with a common 
probability distribution P. The model is assumed to be ferromagnetic, > 
0, i.e. P is supported in The family of random variables u = {uJx,x G 
Z} is an ergodic random field on Z with the space of realizations and 
the probability distribution P = P^. It is known that for every bounded 
realization of the random field u and that for P-a.e. unbounded u the random 
spin system with Hamiltonian (|l]) has a unique limit Gibbs measure for 
arbitrary inverse temperature /5 ^, |^ . To simplify our notation we include 
P into the definition of coupling uj^- 

For a fixed realization of couplings the Ising spin configuration a evolves 
in time through spin flips as specified by the flip rates 



Thus in a short time interval dt the spin configuration a changes to the spin 
configuration 0"*^^^ with probability c{x,a,uj)dt and remains unchanged with 
probability 1 — c(a;, a, It is proved in that this rule defines a 

Markov process, denoted here by 



with state space fl. We assume that the t = distribution of cr'^{t) is the 
Gibbs measure i>i_j. Then is stationary in time and reversible. The 

corresponding stochastic semigroup Tt is self-adjoint on the Hilbert space 
'^oj = '^2(^7 du^). Tt is generated by 



= 5^c(x,a,c.)(/(a(^')) - /(a)), fia) eVcH^, (2) 



c(x, a, uj) 
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X -|- g-Ai(CT,a;) ' 




a'^(t) = G Z,t > 0}, 
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as acting on cylindrical functions V C Tiuj- The operator L{ijj) can be ex- 
tended in 7Y(^ to a self-adjoint (unbounded) operator for P-a.e. P, |^ and 
will be denoted by the same symbol. Let 

Ao = sup{(L^,^),||^|| = l,(V',l)=0} 

denote the upper edge of the spectrum of the operator L{ijj) in the subspace 
orthogonal to {1}. Aq is constant almost surely. 

The goal of our paper is to obtain the long-time behavior for the disorder 
averaged time-autocorrelation function 

^W = (KW>^o(0))p(,)), t^oo. (3) 

Here {■)v{uj) is the average over the the process cr'^(t) under a fixed realization 
a;, and (■) is the average over the distribution P of random couplings. We 
are interested in the case of unbounded couplings when 

P{u^ > K)>0 for any K > 0. 

In this case with probability one the operator L has no spectral gap, see for 
example 0, which implies Aq = 0. 
As our main result we state 

Theorem 1. Let P{uJx > K) > for every K > 0, and 

1 < ((coshtu^)'^) < oo. (4) 
Then for large enough t the following estimate holds 

CA <S{t)<cA (5) 

with positive constants Ci,C2 independent on t. Here 

gi{li) = InP ( u;^ > ^In - ) , g2{li) = 21nP ( u;^ > ^ In — 

with ji G (0,1) for suitable constant c, < c < 1. Gj is the Legendre 
transform of gj, 

Gj{t)= mm (tiJ,- gjiiJ,)), t > 0, j = l,2, 
/xe(o,i) 



where the minimum is taken at fij{t). 

Examples, (i) If P{uJx > m) ~ e^'^" for n ^ oo with k > 4, then one has 

k k 
9i{f^,t) = -\nn, /ie(0,l), = — , 

k 

g2{fi,t) = klncfi, fie {0,1), fi2{t) = j, 

and 

C2(l + t)-''<5(t)<Ci(l + t)-t 
with constants Ci,C2 independent on t. 
(ii) If P{ujx > n) ~ e"""" for u oo with a > 1, then one has 

f/i(/x,t) = -Q^ (in^) , /ie(0,l), 

/ii(t) = af-J (1 + 0(1)), t^oo, 

^72(/i,t) = -2Q^ (^Ini-^ , /iG (0,1), 

/i2(t) = 2a(-j ^ ; (1 + 0(1)), t^oo. 



and 



with some positive constants Ci,C2. 

Our analysis estimates the integrated density of states of the generator 
L in the one-spin sector near zero. Using techniques from the oscillation 
theorem (see [§, |^, for example) we establish a relation between realizations 
of the random couplings and the integrated density of states. This approach 
is exploited already in 0, |[ in the case of bounded u. There the sub-leading 
correction to the exponential decay of S{t) is determined through an analysis 
of the asymptotics of the density of states. For unbounded couplings, how- 
ever, a new mechanism appears resulting in a novel behavior for the spectral 
characteristics of the generator. For bounded couplings the spectrum near 
the upper edge comes from low probability, atypical random couplings, for 
which there are long stretches close to the maximum. They result in a Lifs- 
chitz tail in the integrated density of states. This type of spectrum boundary 
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is called fluctuation boundary. In contrast, as follows from the arguments 
given below, for unbounded couplings the main contribution to the spec- 
trum close to zero comes from rapid oscillations of the couplings over short 
intervals. This behavior of the spectral characteristics of the generator deter- 
mines the leading decay of the disorder averaged auto-correlation function. 
In particular, this implies that the sub-leading decay for bounded couplings 
is unrelated to the leading decay for unbounded couplings. 



2 Reducing subspace, proof of Theorem 1 

The auto-correlation can be rewritten as follows 

Sit) = ((e*^(-)cro,ao)) = ( (e'^^^^Vo, ^o) ). (6) 

We must explain the meaning of the operator Li = Li{uj). Since the work 
of R. Glauber |]10[ it is known that the linear span of "one-point configu- 
rations" {ax, X G Z} forms an invariant subspace for the generator L{uj) of 
(]^). Moreover, the same invariant subspace Hi{uj) C H^j is spanned by the 
functions 

Vx{<J, uj) = coshtUi: ■ ax — sinhux ■ ax-i, x G Z, (7) 

see 0. The functions (0) form the orthonormal basis in Tiiluj). We denote by 
Li(co') the restriction of the generator L{uj) to the invariant subspace Ti.i{uj). 

The operator Li has the following symmetric representation in the basis 
{vx,x G Z} 

LlVx = Ax^x-l^x-l + ^x,x'^^x + ^x,x+lVx+l, (8) 

with 

A -A - ^^\/(^ _, , 

^x,x-l — ^x-l,x — 7Z ^ — 9 ^ ) flz — ISHSlUx > U, 

(1 -aX-i) 



A — —1 



(1 - alal_^) (1 - alal^^] 
We consider new random variables 
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Then 

Using the representation (^) one can introduce the integrated density of states 
iV(Li, dX) for the random operator Li by the truncated operators L^i\r G N, 
defined on a finite-dimensional space of functions Vr of the form 

I, x=—r J 

Let Py,, be the projection on Vr. Then the truncated operator L^^^ = 
Py^LiPy^ is given by the same formula as when x = — r, . . . ,r. We 
denote by > ^ > > • • • > Ag^+i the eigenvalues of the truncated 
operator L^^^ in decreasing order and by k{Lf \ A) the number of eigenvalues 
of L^i^ exceeding A G M. Then from results in 0] it follows that there exists 
a non-random positive measure N{Li,dX) on M, such that with probability 
one 

hm -^kiLP,X) = N{L,,X) 
in the sense of weak convergence of measures, where 

iV(Li,A) = iV(Li,(A,+oo)). 

In addition 

iV(Li,A) = {{ElAK+^)vo,Vo)), (10) 

where {EL-^{dX)} is the spectral resolution of the operator Li. The represen- 
tations (P) to (P) imply (see Lemma 1 below) that the measure N{Li,dX) is 
concentrated on M_, so that N{Li, A) = N{Li, (A, 0)) for negative A. 

Main Lemma. Let A < with \X\ sufficiently small. Then 



N{Li,X) > C2 



P I LJr > - In 



2 clAI 



Cse^^d^'), (11) 



N{L,, A) < C^P (^LJx > i In = Cie^^d^l) (12) 
with positive constants Cj, j = 1,2 and a constant c, < c < 1. 
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The proof of the main lemma will be given in Sections 3 and 4 below. We 
first derive the asymptotic formula (|^) based on the estimates {^), (0). 

Proof of Theorem 1: 

l)The upper bound. Since, see [0, 0, 

a^, = ^D^^y{uj)vy, (13) 

y<x 

with 

Dx,y{uj) = (1 — tanh^ tanhcUy+i . . .tanhu;^., y < x, 

D,,,{cu) = (l-tanh'w,)i/2, 
we have 

x<Q y<0 

< EE «o^S,o>'''((e*"^-., (15) 

x<0 y<0 

2 



(14) 



.x<0 



The representation ( |14| ) together with the condition (^ on the distribution 
of the random variables imply that for any x < 

(Dlo) = ((1 - tanh^o;^)^) (tanh^cu^+i) ■ ■ ■ (tanh^L^o) < k'"^', 

with some < k < 1, so that 

E«o>'^'<C = C(/.). (16) 



Furthermore, for every x, y 

{e'^^'v^.Vyf < {e'^^'v^^v,) {e'^-'Vy^Vy) < {e'^^'v^^v,) . (17) 
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Finally from (|T^ - ([T7|), (]TI1|), and (JT^ we conclude that for large t 
S{t) = {{e'^^ao,ao))<C'{{e'^^v„vo))'^' 



where gi, Gi, and fii(t) are defined in Theorem 1. 

2) The lower bound. By (|^) and we obtain in analogy with the above 
reasoning that 

( (e*'^^fo, I'd) ) = ( (e*^^ ((To cosh cjq — 0"_i sinh cuo), ctq cosh a;o — (T-i sinh cuq)) 
< (cosh^cuo (e*^Vo,ao)) + (sinh^ cjq (e*^'(X-i, (X-i) > 

(18) 

1/2 / 2 \ "^/^ 

)|> < 2(cosh^u;o> ' ((e*^Vo,cro)^^ 



+ 2(cosh2.;osinh2^o>'^'((e*^Vo,a_i)')'^' < fci ((e^^^^o, ^0))'^' 
with some constant ki, where we used the estimate 

(e*^Vo,a_i)' < (e*^Vo,ao) (e*^V_i,a_i) < (e*^Vo,(To). 

Now from (^), (|T^), (0), and (|TT]) we derive for large t the lower bound on 
S{t) as 

= A; / e*^A^(Li,rfA) > ^2 t / e~'^'+<^^^''^d^i 



>C2 



where 5^2? ^52; and /i2(i) are defined in Theorem 1. This completes the proof 
of the theorem. □ 
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3 The estimate of N{Li,X) from below 

Let us fix the configuration uj = {ux, x G Z}. The truncated operator L^i\uj) 
defined above by (§) is given by a Jacobi symmetric matrix of the order 2r + 1 
with positive entries A^^x-i, x = —r + 1, . . . , r. Consequently, for any r the 
operator L^{^ has only real eigenvalues and we can exploit the technique of 
the oscillation theorem in the spectral analysis for L^[\ 

Lemma 1. For every r and f & Vr one has 

0< (-LfV,/) <2l|/f. 
Proof: The proof easily follows from the obvious inequalities 

2^/^2/(1 - Cy)fyfy^l < {l - C y) f^^^ + C y , 
2^Cy{l-Cy)fyfy,, > - {l - C y) - C y . □ 

(r) 

Lemma 1 implies that the operators L\ have only negative real eigenval- 
ues = 1, . . . , 2r + L First we evaluate the function k{L^{\ A) from 
below for A < 0. 

Definition. We call a bond {x,x + 1} regular, if the random variables Cx 
and Cx+i, defined by (||), satisfy the condition 

1 + Cx - Cx+i < \X\. (19) 

Then the following estimate holds. 
Lemma 2. For given A < 

k{L^;\X)>nr{X), (20) 

where TZr{\) is the number of regular pairs, arranged on the interval [— r, r] 
without overlapping. 

Proof: To calculate the number of eigenvalues of L^^'' exceeding A < we 
will exploit the oscillation theorem to the operator — L^^^ and estimate the 
number k{—Li \ |A|) of eigenvalues of —L^^'' not exceeding |A|: k{L^{\X) = 
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|A|). Let {/z(A)} is an eigenfunction of —L^ corresponding to an 
eigenvalue |A|. We define the standard phase V5x('^) by 

, / \ . /x+i(A) 

ctg v?x+i(c^) = ctg v?x+i = r /^x . X = -r, . . . , r - 1. 

Then 

i + a-a+i-|A| - a) i 

ctg 09^,. . 1 = — , , ■ . (21) 

Va+i(i-a+i) va+i(i-a+i) ctg^. 

By the oscillation theorem k{—L'{\ |A|) = m.r(J(A)) + 1, where J(A) < |A| 
is the maximal eigenvalue of —L^ not exceeding |A|, and mr{uj,J{X)) is 
the number of sign changes in the sequence of coordinates {/z(</(A))}, x = 
— r, . . . ,r of the corresponding eigenfunction. Thus m^{J{X)) equals the 
number of sites x G [— r, r] with ctg ip^ < 0, 

mr{J{\)) = #{x G [-r,r] : ctg < 0}. 

Let us consider a regular bond {x,x + 1}. If ctg < 0, then we already 
have a contribution to mr{J{X)) from that bond. If ctg ip^. > 0, then ( pi]) 
and (|19D imply that ctg v^x+i < 0. So in any case we have a contribution to 
mr{J{X)) from each regular bond. Lemma 2 is proved. □ 

Finally by averaging the inequality (^) over realizations u and taking the 
limit r — > oo we have for A < 0, 

iV(Li, A) = hm > 6P(1 + Co - < |A|) (22) 

r-+oo 2r + 1 

with some constant b. We estimate the probability P(l + Co — Ci < |A|) 
under sufficiently small |A| in terms of the distribution P of uj^- 

Lemma 3. For all sufficiently small |A| 



P(l + Co-Ci < |A|) >po 
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(23) 



with constants < po < 1 o,nd < c < 1 . 

Proof: Let us fix some constant h,0 < h < 1, and we denote by 
Po = P{0 < tanhtu^ < h), < po < 1. 
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Then using the representation (j^) for we have for small enough |A| 
P(Co + 1 - Ci < |A|) > P(Co < |A|/2; 1 - Ci < |A|/2) 



1 — a^a?L^ 



< |A|/2; 



1 — a 



1 — a^af 



12 < 1^1/2 



> P(a_i > 1 - Co|A|; < Oq < /i; cii > 1 - co|A|) 



= po[^K>l-£o|A|)]' 
with some cq and < c < 1. □ 



P I > - \n ^ 



" 2 c|A|;_ 



The estimate (|TT]) on N{Li, A) from below follows from (12^) and 



4 The estimate of N{Li,X) from above 

For given A < we denote by 7(A) = iln|i|. Then for any configuration 
u = {ux, X G Z} we consider a decomposition of Z into two sets, 

Z = A^^x U 5^,A 

with 

A^,x = {x G Z : ^, > 7(A)}, 5,,A = {x G Z : < 7(A)}. (24) 
For any r G N we denote by 

Br,uj,x = {x e [-r,r] : max{u;^, < 7(A)} (25) 

5r^,ACfi^,An[-r,r], 

and by Wr^x C K the linear span of functions {vx, x G -6°^ a}- Then the 
operators 



(26) 
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are truncations of Li on subspaces Vr and Wr^x respectively. Since Li is a 
self-adjoint bounded operator, we have by the minimax principle 

kiL^^^\\) > kiLt\\) - #{x G [-r,r],x^ 5,%}, (27) 

where, as above, k{A, A) denotes the number of eigenvalues of the operator 
A exceeding A. 

Lemma 4. For any sufficiently small \X\, A < 0, and for every uj we have 

fc(LS^''^^(cu),A) = 0, (28) 

where L^^'''^\uj) is defined in (^2d^). 

Proof: We consider the bounded configuration 

^ = {^x < e z}, 

coinciding with the configuration uj on B^^^^, 



OJcr. 



Let LKlj) be an operator in Hi^u) given by (H), (^ and corresponding to 
the configuration uj. Our constructions (p4|) - (l25| ) imply that the operator 
(u) is the same as the truncation of the operator LKu) on the same 
subspace Wr,\. As follows from results of [Q] in the case of bounded couplings, 
under the assumption < 7(A) the upper spectrum edge of the operator 
LKu) equals to 

Ao = -1 +tanh27(A) 

for a.e.-configuration cD, so that Aq < |A for small enough A < 0. This esti- 
mate is valid also for any truncation of the operator Lj{uj). Thus no eigen- 
value of the operator Lj' ^ (uj) or (uj) can be greater than A. □ 



By (|27|) and (^8]) we have the following estimate 

k{Lt\u),X) < #{^e[-r,r],a;^i?,%} 

< C#{xe [-r,r],xe A^,a} (29) 

with some constant C. Applying, as before, the ergodic theorem to the 
inequality (0), we obtain the estimate ([T2|) on N{Li, A) from above, 

NiL„ A) = hm iM^lM < c^P U > i In ^ 
r^oo 2r -|- 1 V 4 A 



12 



References 

[1] R. A. Minlos, A. G. Trishch, The complete spectral decomposition of a 
generator of Glauber dynamics for one-dimensional Ising model, Uspechi 
Mathem. Nauk 49, 209-210 (1994). 

[2] S. Albeverio, R. Minlos, E Scacciatelli, E. Zhizhina, Spectral analysis of 
the disordered stochastic 1-D Ising model, Commun. Math. Phys. 204, 
651-668 (1999). 

[3] B. Zegarlinski, Strong decay to equilibrium in one- dimensional random 
spin systems, J. Stat. Phys. 77, 717-732 (1994). 

[4] G. Gielis, C. Maes, The uniqueness regime of Gibbs fields with un- 
bounded disorder, J. Stat. Phys. 81, 829-835 (1995). 

[5] T. Liggett, Interacting Particle Systems, Springer- Verlag, Berlin, 1985. 

[6] S. A. Gredeskul, L. A. Pastur, Behavior of the density of states in the 
one-dimensional disordered systems near the spectrum bounds, Teoret. 
and Matemat. Physika 23, 132-139 (1975). 

[7] E. Zhizhina, The Lifshitz tail and relaxation to equilibrium in the one- 
dimensional disordered Ising model, J. Stat. Phys. 98, 701-721 (2000). 

[8] E. Zhizhina, Spectral analysis of an one-dimensional stochastic Ising 
model with random potential: asymptotics of the time auto-correlation 
function. Trans, of Moscow Math. Society 64, (2002), to appear. 

[9] L. Pastur, A. Figotin, Spectra of Random and Almost-Periodic Opera- 
tors, Springer- Verlag, Berlin, 1991. 

[10] R. Glauber, Time dependent statistics of the Ising model, J. Math. Phys. 
4, 294-307 (1963). 



13 



